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Abstract 

A  cracked  metallic  coupon,  called  crack  gage,  is  being  considered 
as  a  device  for  monitoring  crack  growth  in  aircraft  structures.  For 
this  purpose,  a  stress  intensity  factor  solution  for  the  gage  has  to 
be  known. 

This  study  provides  stress  intensity  factor  solutions  for  two  basic 
geometric  configurations  subjected  to  prescribed  displacements. 

1.  Edge  cracked,  trapezoidal  shaped  gages  of  uniform  thickness. 

2.  Center  cracked  gages  with  varying  or  stepped  thickness. 

For  the  trapezoid,  the  influences  c£  changing  the  length  of  the 

cracked  edge,  while  other  edge  remains  constant,  and  vice  versa,  were 
investigated.  The  results  obtained  do  not  show  significant  beneficial 
changes  in  stress  intensity  factor  fo?;  the  range  of  parameters  con¬ 
sidered  over  those  of  rectangular  gages . 

Stress  intensity  factors  were  ^determined  for  stepped  gages  of 
various  geometries.  Various  thickness  ratios,  length  ratios  and  aspect 
ratios  were  considered,  including  the  specific  geometries  of  two  gages 
now  under  development.  In  each  case,  the  stress  intensity  factor  was 
determined  as  a  function  of  crack  length. 

Stress  intensity  factors  were  found  to  increase  as  the  cracked 
center  portion  was  made  thinner,  and  as  the  length  of  the  outer  section 
was  increased.  The  stepped  gage  was  found  to  have  the  potential  for 
tailoring  the  stress  intensity  factor.  Finally,  analytical  results 
of  stress  intensity  factor  were  found  to  show  good  agreement  with  exper¬ 


imental  results. 


EFFECT  OF  GEOMETRIC  VARIABLES  ON  STRESS 


INTENSITY  FACTORS  FOR  CRACK  GAGES 

I .  Introduction 

Background 

The  concept  and  the  requirements  for  safety,  durability  and  man¬ 
agement  of  military  aircraft  structures  are  summa  -ized  by  Coffin  and 
Tiffany  (Ref  1) .  As  reflected  there,  the  durability  analysis  as  well 
as  scheduled  maintenance  for  the  structure  are  based  on  some  assumed 
usage  of  the  fleet.  Usually  the  actual  usage  of  an  individual  air¬ 
craft  and  the  fleet  will  differ  from  the  usage  for  which  the  aircraft 
was  designed.  In  order  to  provide  updated  usage  data  for  force  manage¬ 
ment,  MIL  STD  1530A(11)  (Ref  2:23-25)  requires  establishing  aircraft 
structures  tracking  programs .  Such  a  program  provides  baseline  oper¬ 
ational  load  spectra  for  average  fleet  usage  as  well  as  an  Individual 
Airplane  Tracking  (IAT)  program,  which  provides  crack  growth  information 
concerning  each  individual  aircraft .  The  IAT  is  the  feedback  link 
necessary  for  the  scheduling  of  structural  maintenance . 

A  typical,  current  technology,  IAT  program  is  that  of  the  F-4 
aircraft.  This  program  was  developed  as  part  of  the  Damage  Tolerance 
Assessment  study  for  F/RF-4  C/D  and  F-4E(S)  aircraft  (Ref  3,  4) .  The 
program  uses  as  tracking  devices  counting  accelerometers  installed  on 
each  aircraft  and  flight  data  recorders  installed  on  12%  of  the  fleet. 

A  crack  growth  analysis  is  performed  for  critical  structural  locations. 
Stress  exceedances  spectra  for  each  critical  location  and  for  different 
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usage  categories  were  developed  from  the  flight  recorded  data.  These 
spectra  along  with  counting  accelerometer  exceedances  readings  are 
used  to  determine  the  stress  history  at  critical  locations .  Using 
these,  a  crack  growth  analysis  for  an  individual  aircraft  may  be  per¬ 
formed  . 

This  tracking  system,  described  in  more  detail  in  other  places 
(Ref  3,  4,  5,  6:64-66),  is  very  indirect  in  defining  stresses  and  in¬ 
volves  many  sources  for  error.  Further,  the  enormous  amount  of  data 
to  be  collected  in  order  to  get  the  stresses  adds  to  program  cost.  A 
device,  the  crack  gage,  has  been  proposed  which  might  provide  a  more 
direct  and  efficient  technique  for  use  in  IAT  programs .  It  has  been 
estimated  that  using  crack  gages  in  the  F-4  XAT  program  might  save 
7.7  million  dollars  in  five  years  (Ref  6:79). 

Crack  Gage  Concept 

A  preliminary  concept  of  a  crack  gage  was  patented  by  Smith  (Ref  7) 
in  1976 .  A  more  practical  concept  for  using  a  crack  gage  for  monitoring 
structure  crack  was  introduced  by  Crane,  Gallagher  and  Grant  (Ref  8) . 
Figure  1  shows  the  concept  of  the  crack  yage .  It  consists  of  a  pre¬ 
cracked  thin  metallic  sheet  with  ends  bonded  to  a  structure .  Since  the 
structure  and  the  gage  will  be  exposed  to  the  same  load  history,  it  is 
possible  to  relate  the  crack  growth  in  the  structure  to  that  in  the 
gage.  The  readily  observed  crack  growth  in  the  gage  can  then  be  used 
as  a  prediction  of  crack  growth  in  the  structure. 

Efforts  in  various  directions  were  made  in  evaluating  the  crack 
gage  concept  since  it  has  been  introduced.  Torvik  (Ref  9,  10)  developed 
an  analytical  solution  for  the  stress  intensity  factor  for  displacement 
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controlled,  edge  cracked  gages.  Ori  (Ref  11)  and  Ashbaugh  and  Grandt 
(Ref  12)  evaluated  the  concept  using  both  edge  cracked  and  center 
cracked  gages . 

Stepped  center  cracked  gages  are  now  being  investigated  under 
two  U.S.  Air  Force  contracts,  one  with  Boeing  VJichita  (Ref  13)  and  one 
with  McDonnell  Douglas  Corporation  (Ref  14) .  The  contract  with  McDonnell 
Douglas  will  provide  evaluation  of  a  center  cracked  stepped  gage  on  a 
F-4  aircraft  during  a  full  scale  fatigue  test. 

Statement  and  Scope  of  Problem 

It  has  been  observed  (Ref  9,  10)  that  relatively  high  stress  inten¬ 
sity  factors  for  an  edge  cracked  rectangular  gage  may  be  obtained  if 
the  gage  is  long,  but  the  magnitude  then  varies  considerably  as  the 
crack  grows,  on  the  other  hand,  constant  stress  intensity  factor  vs 
crack  length  may  be  obtained  for  relatively  short  gages,  but  the  gage 
sensitivity  is  then  very  lov. .  For  analytical  purposes ,  such  as  pre¬ 
dicting  the  remaining  life  of  a  structure,  a  constant  stress  intensity 
factor  vs  crack  length  is  desired.  It  was  suggested  that  an  edge-cracked 
gage  of  a  trapezoidal  planform  might  satisfy  the  requirement  for  a 
higher,  but  constant  Kj  vs.  crack  length,  since  the  "effective"  length 
of  the  gage  could  vary  as  the  crack  grows.  The  stepped  gage  was  also 
considered  since  the  thicker  section  serves  as  an  amplifier  of  the 
strain  in  the  cracked  section  from  the  far  field  value .  The  gain  may 
be  controlled  by  varying  the  thickness  ratios  and  length  ratios  of 
the  segments  of  the  gage. 

The  gages  that  had  been  developed  for  the  activities  previously 
mentioned  were  basically  developed  empirically.  It  became  apparent 


that  an  analytical  method  that  could  be  used  to  obtain  stress  intensity 
factor  solution  for  stepped  center  cracked  gages  was  required.  Gages 
with  different  stress  intensity  factor  sensitivities  are  required  so 
that  the  gage  can  be  matched  to  the  specific  usage,  i.e.,  applying 
gages  of  very  high  sensitivity  to  structures  to  be  subjected  to  a  large 
number  of  cycles . 

In  this  study,  a  means  of  obtaining  analytical  predictions  of  the 
stress  intensity  factor  for  the  two  geometries,  trapezoidal  planform 
and  stepped  gages,  is  developed.  These  methods  are  applied  to  the  design 
of  crack  gage  and  the  results  for  the  stepped  gage  were  experimentally 
verified . 
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II.  Stress  Intensity  Factor  Solution  for  a  Cracked 
Flat  Sheet  Subjected  to  Mixed  Boundary  Conditions 

Introduction 

A  method  utilizing  Williams'  (Ref  15,  16)  stress  functions  along 
with  energy  principles  was  developed  by  Torvik  (Ref  9,  10)  to  obtain 
a  stress  intensity  factor  solution  for  a  finite,  cracked  elastic  sheet. 
Solutions  were  provided  for  a  rectangular  edge  cracked  sheet  subjected 
to  mixed  boundary  conditions.  Part  of  the  boundary  was  subjected  to 
prescribed  tractions  as  shown  in  Figure  2A,  herein  referred  to  as  the 
type  1-mixed  boundary  condition.  The  solution  for  a  center  cracked 
sheet,  which  is  the  center  section  of  the  stepped  cracked  gage,  led  to 
a  different  mixed  boundary  condition,  herein  referred  to  as  type  2, 
and  shown  in  Figure  2B.  For  this  problem,  a  part  of  the  region  of  the 
physical  surface  is  subjected  simu) taniously  to  prescribed  tractions 
as  well  as  to  displacements.  A  brief  outline  of  the  method  developed 
by  Torvik  (Ref  9,  10),  for  the  type  1-mixed  boundary  condition  problem, 
and  the  necessary  extension  for  the  type  2-boundary  condition  problem, 
will  be  presented  here . 

Method  of  Solution 

Consider  a  thin  elastic  sheet  for  which  plane  stress  conditions 
hold,  with  a  crack,  starting  at  an  edge.  A  coordinate  system  is  set 
with  origin  placed  at  crack  tip.  The  crack  faces  are  located  at  9=±tt 
OSrSa,  and  are  free  of  tractions.  Williams  (Ref  15,  16)  has  shewn 
that  for  an  infinite  domain,  stress  fields  can  be  obtained  by  using  a 
stress  function  of  the  form 

X(r,  0,  X)  =  rX+1  F(G,  X)  (1) 


G 


A  =  m/2  for  m  =  1 ,  2 ,  3 ,  . . . 


(2) 


The  stress  function  may  be  divided  into  even  and  odd  parts .  since 
the  geometries  of  interest  for  this  study  are  symmetric  in  respect  to 
X  axis,  only  the  even  part  of  the  stress  function  is  of  interest.  In 
addition,  displacements  corresponding  to  the  stress  function  may  be 
found  from 
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2yU0  =  -  i  +  (l  -  S)r2  ^ 
r  30  8r 

jj  =  shear  modulus  s  =  V/l+V 


(4) 


V  =  Poisson  ratio 


2  3  ,  S'K 

v  X  =  —  (r-^) 
3r  39 


(5) 


The  even  parts  of  x  2111(1  ^  >  as  well  as  the  corresponding 

stresses  and  displacements,  are  shown  in  Appendix  A.  Each  stress  and 
displacement  is  made  up  from  an  even  and  an  odd  part.  As  stated  by 
Torvik  (Ref  9:3-4,  10),  the  stresses  o  and  the  displacements  U  may 
be  represented  in  the  form: 
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where  each  am  satisfies  field  equations  and  the  traction  free  require' 
ment  on  the  crack  face,  and 
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Each  of  Um  satisfies  strain  displacement  equations  and  stress 
strain  laws,  and  the  arbitrary  coefficients  dg,  cure  related  to  the 
stress  and  displacement  functions  coefficients  presented  in  Appendix 
A  as  Am  through  the  relation 

<*1  ■  Ai,  d2  *  -A2,  d3  ■  -A3>  d4  ■  A4,  d5  =  A5,  d6  =  -Ag,  etc. 

In  order  to  find  a  stress  intensity  factor,  defined  as, 

Kj  m  lim{/27Tr '  OqqCG  =  0,r)}  (8) 

r-K) 

it  will  be  required  to  find  the  stress  Oqq  for  a  certain  boundary 
condition  problem.  More  specifically,  the  coefficients  ^  in  Eq  (6) 
are  to  be  obtained. 

After  applying  the  limit  procedure,  the  only  term  to  remain  is 
the  one  multiplied  by  d^  .  All  other  terms  vanish  since  they  include 
i ?  expressions ,  pt  1/2 .  Hence , 

K_j.  -  /zP  Ax  =  -/IP  d1  (9) 

The  solution  of  the  coefficients  for  a  mixed  boundary  condition  of 
type  1  was  given  by  Torvik  (Ref  9;20-21,  10)  utilizing  Reissner's  prin¬ 
ciple.  The  interpretation  of  Reissner's  principle  by  Fung  (Ref  17:299- 
300)  is  that  the  surface  S  enclosing  the  domain  V  ,  can  be  divided 
into  Sip  and  Su  ,  where  tractions  are  prescribed  on  ST  and  displace¬ 
ments  on  Su  .  A  more  general  interpretation  will  be  presented  here. 

In  the  type  2  mixed  boundary  condition,  a  surface  region  on 

the  surface  S  can  be  subjected  simultaniously  to  prescribed  tractions 
and  displacements  as  well.  Hence,  it  may  be  considered  a  subset  of  Sfj. 


as  well  as  S„  ,  with  the  relations : 


s  -  ST  +  su  (10) 

=  +  ^Tu  (11) 

su  “  Cu  +  (12) 

where  Gj.  is  the  region  where  only  tractions  are  prescribed,  and  cu 
is  the  region  where  only  displacements  are  prescribed. 

Consider  the  functional  (essentially  Reissner's  function) 

J  m  Jl  W(e±j)  -  F^jdv  -  fa ijfe^  -  1/2 (U± j  +  Uj#iJ]dv  - 

v  v 

-  I  Ti*  u^ds  -  I  cr± j  Vj  (Hi  -  u±*)ds  (13) 

^ST  \ 

Where  stresses ,  displacements  and  strains  are  to  be  independent  vari¬ 
ables,  and  =<J ji .  The  elastic  equilibrium  will  be  obtained  by  setting 
J  “  0,  or  j  will  be  stationary  with  J  =  0.  The  a_.  are  introduced  as 


Lagrange  multipliers. 

6j  =  0  =m|5-  6e.  . 

J  3eij  13 


3w  f 

-  piSui>*'  -  Jto  ijlVj  -  1/2(0,  .,  +  o3ii))a»  - 

>  y  V 

-/  -  l/2(6uif j  +  6Ujfi)]dv  -j  T±*6u^d s  - 

v  ^  St 

-  j  6(vj0ij)  ‘(^  -  U^*)  dS  -  j  VjaijlS  Uids 


Rearranging  and  applying  the  divergence  theorem  to  third  integral,  we  find 

5j  -  °  mJ( f£  -  0ij)6eijdv  -f(F±  +  ^-)<5u.dv  -  f 6a  ^e..  - 
v  “V  ■*  '  v 

-  l/2(U1#j  +  Uj>A)]dv  +  fvj  OijfiUidS  -  j  T^fiU^dS  - 

S  ST 

-  j  (0i  -  U^JiSTjdS  -  J  T^l^dS  (15) 
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Realizing  that  the  first  three  integrals  represent  the  stress-strain 
law,  equation  of  equilibrium  and  strain  displacement  equation  respec¬ 
tively,  they  may  be  omitted  from  Eq  (15)  if  we  choose  the  stresses  to 


be  represented  by  Eq  (6) ,  and  displacements  by  Eq  (7) ,  along  with  the 
corresponding  strain. 

Eq  (15)  reduces  to: 

f^o^sv.as  ■L  T^fiUjdS  - j  (U±  -  Ui*)6TidS  -J  T^l^dS-0  (16) 

/  (Ti  -  T±*) 6UidS  -  /  (Ui  -  Ui*) 6TidS  =  0  (17) 

'  S„  J  s 

T  U 

Applying  Eqs  (11)  and  (12) 

/  (T±  -  Ti*)6uidS  -  (U  -  Ui*)*6Tids  + 

•'/Cm  Jc  1 

T  u 

+1  [(T(0-  T(U)-S0(i)  '  ID(1)  '  u(j*’-ST(j,]as  -  0  (18> 

Td 

where  j^i,  and  the  parentheses  indicate  suspension  of  the  summation 
convention.  On  Cpy  ,  one  component  each  of  displacement  and  traction 
are  given. 

Eq  (18)  may  be  used  to  find  a  Kj  solution  for  problems  with  mixed 
boundary  condition  of  type  2. 

Use  now  the  stresses  and  displacements  as  defined  in  Eq  (6)  and  (7) 
to  find  variations  of  stress  and  displacement. 


5oij  - 


f^d-a  « 

P  ^ 


11 


2y 


6u. 


6u„ 


M 


6dp 


p=l 


U  P 
r 


(20) 


Substituting  Eqs  (6) ,  (7) ,  (19)  and  (20)  into  Eq  (18) : 


M 


f  ttVj  X  Vijm  -  lyE  6dPU iP]dS  ‘ 
SFi  P=1 

■f  l(^.dmuim  -  2yu±*)v.y  6dpa..P]ds-+ 

'  cu  n^L  P^l 


M  M 

f  u^L  VwrWbt-',1*0™' 

J C  3  m=l  ^  P=1 

Tu 


M  A 

-  3T  'I,  V(j)  -  2u°tj7vi^i  Mp"i  $as  ‘  ° 


(21) 


Since  6 dp  is  arbitrary,  make  the  choice  6dp  =  0  ,  p?*q,  6 dp  -1 

if  p=q  and  then: 


M  M 

f  <h £  w>°iq  ^  -/  >  I?»°iB)vjaijqas  ■ 

\hen?i  Jc  ^ 


c  C 

U  Tu 


■J  Ti**UiqdS  -j  [Ti*-UiqdS  -  2yUi*Tiq]dS  -  2)j  U^-TjSdS 


Tu 

Eq  (22)  may  be  represented  as : 
M 


(22) 


E 


m**l 


VS 


(23) 


mq 


where  A^q  is  a  symmetric  matrix. 
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t 


V[  -  Jc  “Vrr"  +  Vr0m)ur  +  'Ve/  +  Ve9)ue]ds  • 

T+  Tu 


Di0-(vA;  +  v0c00q)  +  or“(vrarrH  +  v^)] 

CU+CTu 


ds 


(24) 


The  reader  is  reminded  that  only  one  product  of  force  and  displacement 
is  to  be  evaluated  in  the  contribution  of  to  each  of  these  integrals . 

Fq  =  J  (T*rUrq  +  TQ*Ueq)dS  -  2\f  (U0*Teq  +  Ur*Trq)dS  _ 

■'C. 


-  j  [Tr*Urq  +  Te*U0q  -  2y(U*TQq  +  U*Trq)]dS 

Ctu 


(25) 


Since  the  displacements  are  prescribed  elsewhere ,  the  crack  tip 
has  to  be  allowed  to  translate.  Thus,  as  stated  (Ref  9:20,  10),  we  may 
introduce  a  translation  term  into  the  expansion,  which  will  be 


2yur 

0 

*  Vi  - 

dgcos© 

(26) 

2yu0 

J  0 

"  Vi  - 

-dQsinO 

(27) 

with  orr°  ,  a00°  , 

ar0°  being 

zero. 

By  solving  Eq  (23)  along  with  Eqs  (24)  and  (25)  for  the  coefficients, 
a  solution  for  a  finite,  two  dimensional,  cracked  elastic  sheet  with 
mixed  boundary  conditions  of  type  2  is  obtained.  Note,  that  by  elimi¬ 
nating  the  integral  on  region,  we  get  a  type  1  mixed  boundary 

problem,  and  A^g  and  Fg  will  have  the  same  form  as  shown  (Ref  9:23,  10) 
previously. 

Since  the  particular  gage  geometries  for  this  study  are  the  trap¬ 
ezoidal  gage  and  the  stepped  gage ,  the  formulation  for  getting 
and  Fg  for  these  geometries  will  be  shown. 
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"><»er 


The  basic  boundary  conditions  for  the  trapezoidal  gage  arc  essen¬ 
tially  as  shown  in  Figure  2a,  and  the  formulation  was  presented  (Ref  9: 

23,  10) . 


=  2y 


[  (UX*COS0 


U 

uy  ‘ 


inC)  (vrarjr  j 


+  v  a 


Cr 


l)  + 


+ 


(-Ux*sin0  +  U  cosG) • (vrOrp  + 


V„G 


9  66 


L)  ]aS 


(28) 


When  evaluating  for  the  discussed  geometry,  crack  faces  are  free 

of  traction.  The  first  integral  of  Eq  (24)  will  be  evaluated  on  edged 

1  and  3,  while  second  integral  will  be  evaluated  on  edge  2. 

The  case  with  boundary  conditions  shown  in  Figure  2b  represents  the 

center  portion  of  the  stepped  gage .  Edges  1  and  2  are  part  of  CT  , 

where  tractions  are  prescribed,  while  edge  3  is  part  of  CTu  . 

Taking  into  account  boundary  conditions,  only  a  part  of  the  integral 

CTu  exists  on  edge  3.  Pus  to  boundary  condition  presentation,  A^g 

and  F  will  be  represented  in  cartesian  coordinates,  and  Eq  (24) 

Si 

becomes 


^mq 


-I 

'Cm 


Ty'uJjdS 


f 


-Tu 


UIt>Tq]dS 
x  x 


(29) 


With  Cu  =  0  ,  and  with  the  prescribed  variables  on  CTu  substituted 
into  Eq  (25) ,  we  get 


F 


q 


+  T  *uq)ds 
y  y 


(30) 


Those  last  two  equations  will  be  used  to  solve  the  stepped  crack  gage 
stress  intensity  problem. 
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III.  Stress  Intensity  factor  Solution  for  A  Trapezoidal , 
Edge  Cracked  Crack  Gage  with  Prescribed  End  Displacements 


The  edge  cracked  rectangular  gage  has  a  tendency  to  give  low  stress 
intensity  factor  values  for  short  gages,  while  the  sensitivity  improves 
for  relatively  long  gages  (Ref  9:34-46;  10).  This  has  also  been  ob¬ 
served  in  experiments  (Ref  6:73).  The  disadvantage  of  long  gages  is 
that  Kj  does  not  remain  constant  and  even  decreases  as  c rack  grows . 

This  property  is  not  desired  for  a  crack  gage.  It  was  thought  that 
giving  the  gage  a  trapezoidal  shape  as  shown  in  Figure  3,  a  better  con¬ 
trol  on  the  Kg  variation  with  crack  growth,  and  a  higher  Kg  might 
be  achieved. 


The  suggested  configuration  consists  cf  a  thin  gage  bonded  to  a 
structure  along  the  line  AE.  The  structure  is  loaded  by  stress  Js 
and  displacements  are  transferred  to  the  gage  through  line  AB.  The 
stress  intensity  factor  solution  for  this  case  was  obtained  by  using 
the  method  snown  in  Section  II.  Using  Eg  (28)  along  with  Eq  (24), 

Eq  (23)  gives  a  solution  for  the  coefficient  dg  . 

The  prescribed  displacements  are 


U*  =  -xv  ^ 
Es 


where  y  =  Ha  -  tga*  (b  -  x)  and  v=.3  .  (33) 

The  matrix  A^q  was  evaluated  by  a  computer  program  (see  Appendix 
B)  using  an  expansion  of  thirty-six  coefficients,  including  these  for 
rigid  body  displacements .  The  integration  was  performed  over  300 
boundary  points.  The  variables  for  this  study  were  chosen  to  be  the 
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Structure 


Cs 


Figure  3.  Trapezoid  Crack  Gage 


angle  a  ,  which  was  assumed  to  control  Kj  changes  as  the  crack  grows, 
and  also  to  control  the  magnitude  of  Kj  ,  since  the  gage  length  dimen¬ 
sions  are  changing.  Former  results  for  a  rectangular  gage  (Ref  9:34, 

10)  indicated  that  a  gage  with  an  aspect  ratio  of  2  Ha/b=2  might  be  a 
good  one  for  the  study,  since  for  that  gage  was  fairly  constant 

as  crack  length,  a  ,  increased. 

Results  and  Discussion 

Non-dimensional  stress  intensity  factor  results  for  Hc/b=l  are 
presented  in  Figure  4.  It  can  be  observed  that  changes  in  a  in  the 
range  of  160°  -  200°  have  a  very  little  influence  on  the  magnitude  of 
Kj  .  The  general  behavior  of  Kj  vs.  a  improves  a  little  and  for 
a  =  200°,  KI/ast'i>  is  almost  constant  on  the  range  .4<a/b<.8. 

Results  for  Ha/b  *  1  are  presented  in  Figure  5.  It  can  be  observed 
that  in  comparison  to  former  results  the  changes  in  Kj  are  relatively 
large  as  a  changes  from  165°  to  195°.  This  can  be  related  to  the 
fact  that  under  the  displacement  loading,  the  JCj.  will  be  more  sensitive 
to  length  changes  on  Hc  than  on  Ha  since  the  former  directly  affects 
the  crack  opening,  and  Kj  may  be  considered  proportional  to  crack 
opening . 
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.2  .3  .4  .5  .6  .7  .8 

Figure  5.  Stress  Intensity  Factor  for  Tiapezoid  Crack 

Gage  Ha/b  =  1. 
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IV.  Stress  Intensity  Factor  Solution  for  a_  Stepped 

Center  Crack  -  Crack  Gage ,  with  Prescribed  End  Displacements 

The  stepped  center  crack  gage  concept  is  described  in  Figure  6. 

The  gage  consists  basically  of  three  segments.  The  center  segment  of 
reduced  thickness  which  includes  the  crack,  and  two  thicker  end  seg¬ 
ments  . 

Part  of  each  thicker  segment  is  bonded  to  the  structure,  and  a 
part  (of  length  L0)  is  free .  It  is  assumed  that  the  gage  is  damped 
to  the  structure  along  lines  AB  and  GH.  This  assumption  is  only  an 
approximation  since  the  actual  load  transfer  from  the  structure  to  the 
gage  should  be  modeled  as  loads  transferred  through  springs,  where  the 
adhesive  layer  is  the  spring  element.  However,  the  exact  model  was 
beyond  the  scope  of  this  study.  Displacements  of  line  AB  relative  to 
line  0-0  are  assumed  to  be  transferred  completely  to  the  gage  with  r.o 
losses  in  the  bond  line,  i.e.,  the  bond  line  does  not  rotate. 

Due  to  the  symmetry,  the  center  segment  of  the  gage  needs  to  be 
solved  only  for  a  quarter  of  the  plate,  as  shown  in  Figure  7.  The 
entire  outer  segment,  shown  in  Figure  8,  was  analyzed  in  order  to  obtain 
better  accuracy  with  the  finite  element  analysis. 

Methods  of  Solution 

Two  analytical  methods  will  te  used  to  obtain  the  complete  solution. 
The  Finite  Element  Method  will  be  applied  to  solve  forces  and  displace¬ 
ments  of  the  outer  segment  of  the  gage  and  the  method  described  in 
Chapter  II  will  be  used  to  solve  for  stress  intensity  factors,  stresses 
and  displacements  for  the  center  segment.  Bowie  et  al  (Uef  18:767-77^) 
introduced  also  the  idea  of  getting  stress  intensity  factors  by  using  a 
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Figure  6 . 


Center  Cracked,  Stepped  Gage  Geometry 


stitching  technique.  They  used  one  analytical  method  to  solve  stresses 
and  displacements  for  cracked  and  uncracked  part.  By  stitching  the 
solutions  stress  intensity  factors  were  found  for  center  cracked  uni¬ 
form  thickness  plate  subjected  to  uniform  stress  loading. 

Since  the  line  CD  is  actually  the  connecting  line  between  the 
two  segments,  we  will  demand  displacements  and  forces  produced  by  the 
two  methods  to  be  equal  on  this  line.  Let  the  line  CD  be  discretized 
into  2n  elements,  each  of  length  Dx.  Each  element  will  be  represented 
by  a  nodal  point,  considered  to  be  at  the  center  of  the  segment,  at 
which  forces  and  displacements  are  assumed  to  be  applied. 

The  outline  for  the  solution  is  as  follows: 

{U*}  =  {Uc}  +  {U0}  (34) 

where  U*  ,  Uc  are  the  prescribed  relative  displacement  of  line  AB 
(Figure  6)  and  relative  displacement  of  line  CD  with  respect  to  inertial 
frame  x,  y,  respectively.  Uq  is  the  relative  displacement  of  line  AB 
with  respect  to  CD.  Using  flexibility  matrix  of  center  segment  [sc]  , 
and  flexibility  matrix  of  outer  segment  [s0 ]  ,  we  have 


{Uc}  - 

tsc] 

{pc> 

(35) 

{UQ>  « 

[Sol 

{p0} 

(36) 

On  line  CD,  the  forces 

{Po} 

=  {pc} 

(37) 

After  substitution 

{U*}  -  [[So] 

+  [sc] 

){pc> 

(38) 

All  the  forces  and  displacements  are  expressed  by  x  and  y  direction 
components  consequently  all  matrices  are  partioned  and  the  number  of 
rows  or  columns  is  twice  the  number  of  nodal  points  under  consideration. 
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Suppose  a  unit  load  is  applied  to  each  element  of  the  center  seg¬ 
ment,  successivley .  Using  the  technique  described  in  Chapter  II,  the 


coefficients  of  the  expansion  may  be  solved  for  each  loading  case. 

Prom  the  coefficients  solved  for  a  vnit  load  at  nodal  point  k,  the 
first  is  of  particular  interest  and  will  give  the  stress  intensity 
factor  for  unit  load  at  nodal  point  k  .  When  we  repeat  the  procedure 
for  all  nodal  points,  we  may  construct  an  array  {d}  whose  k'th  element 
is  the  contribution  of  load  at  point  k  to  the  stress  intensity  factor. 
The  stress  intensity  factor  for  the  center  segment  under  uniform  dis¬ 
cretized  unit  load  is  obtained  by  summing  the  elements  of  {d}.  If  the 
load  at  each  nodal  point  is  different,  the  stress  intensity  factor 
solution  may  be  written  as  s 

KX  -  27r{D}T{Pc}  (39) 

or, 

Kz  =  2n{D}r  {[Sc]+  fc0}]"1fU*}  (40) 

Thus ,  the  solution  for  the  stepped  gage  with  prescribed  end  dis¬ 
placements  requires  that  the  flexibilities  [  Sc]  and  [sQ]  be  determined. 

A  prescribed  stress  T*  ,  on  k'th  element,  Dx  long,  may  be  con¬ 
verted  into  a  "concentrated"  stress  at  the  nodal  point  xk  producing 
a  unit  force  per  unit  thickness  by  defining  a  loading  function. 


T*  -  fi(x  -  xk) 


T*  -  <5(x  -  Xfc) 


(41) 


ao  as  to  produce  an  unit  force  in  each  coordinate  direction  on  any 
element,  centered  on  xk 
x^+px/2 

5(x  -  xk)dx  »  1  (42) 


xk-Dx/2 
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For  simplicity  a  unit  thickness  is  considered.  The  vector  Fg 
can  be  computed  by  Eq  (30) ,  with  displacements  evaluated  at  x^  . 

The  matrix  A^g  can  be  evaluated  by  Eq  (24)  and  the  system  of  equations 
can  be  solved  for  the  coefficients  by  using  Eq  (23) .  With  the 

set  of  coefficients  known,  the  actual  displacements  at  every  other 
nodal  point  1  can  be  evaluated.  By  repeating  the  procedure  for 
k  =  l...n,  we  contruct  the  flexibility  matrix  [3C]  ,  along  with  the 
array  {d!  which  is  constructed  from  the  coefficients  d^  for  each  of 
k  load  points.  The  flexibility  matrix  [Sc]  is  actually  factored  such 
as  to  account  for  actual  thickness  ratio  between  center  and  outer  seg¬ 
ments  of  gage.  The  outer  segment  flexiblity  matrix  is  contructed  by 
using  "ANALYZE"  (Ref  19) ,  a  finite  element  analysis  computer  program 
developed  at  the  Air  Force  Flight  Dynamics  Laboratory.  The  matrix  was 
established  by  constraining  one  end  of  the  segment  completely  to  simu¬ 
late  the  bond  line,  and  the  other  end  was  loaded  with  unit  forces  at 
discrete  points . 

Numerical  Evaluation 

The  problem  introduced  for  the  center  segment  was  solved  ty  com¬ 
puter  programs  which  are  described  in  Appendix  B.  The  objective  of  the 
pzograxr.  was  to  evaluate  the  matrix  A^  ,  given  in  Eq  (29)  ,  to  solve 
Eq  (23)  for  the  coefficients,  and  to  construct  the  flexibility  matrix 
[sc]  in  the  manner  previously  described. 

The  matrix  A^q  was  evaluated  by  a  numerical  integration  along 
edges  1,  2,  3,  using  Simpsons'  third  rule.  The  boundary  was  divided 
into  300-400  points.  A  slight  sensitivity  to  the  number  of  integration 
points  was  observed  (for  short,  Lc)  when  elements  of  different  lengths 
were  used  to  perform  the  integration  on  separate  edges . 
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Hie  number  of  terms,  M,  used  for  the  expansion  was  18.  The  number 
of  nodal  points  used  on  edge  2  was  20.  Consequently,  the  flexibility 
matrixes  were  dimensioned  40x40,  partioned  into  x  and  y  components. 

The  flexibility  matrix  of  the  outer  segment  was  developed  for  a 
whole  outer  plate,  and  the  corresponding  40x40  matrix,  J^SQ  ,  was 
extracted . 

The  prescribed  displacement  U*  was  chosen  to  be  1  in  longitudal 
y  direction  and  zero  in  x  direction.  As  has  been  shown,  releasing 
the  transverse  constraints  would  have  little  effect  on  results  (Ref  9: 
32,  10) ,  since  the  stress  intensity  factor  is  driven  primarily  by 
displacments  perpendicular  to  crack  plane - 

Stress  intensity  factors  were  computed  for  center  segments  having 
lengths : 

Lc  =--  .4;  .6;  .8;  1.0;  1.25 

The  outer  segments  having  lengths 

L0  =  1.0;  1,267;  1.6;  2.0;  4.0 

and  for  the  thickness  ratios 

tc/t0  =  .2;  .5;  .75 

The  width  in  all  cases  was  taken  to  be  b=2.0. 

Stress  intensity  facors  were  not  calculated  for  all  combinations  of 
variables . 

Results  anjd  Discussion 

A  sample  of  the  results  of  a  convergence  study  for  the  stress  in¬ 
tensity  factor  solution  of  the  center  segment  with  12,  18  and  24  terms 
is  presented  in  Table  I.  No  significant  changes  for  more  than  18 
terms  were  observed.  The  number  of  nodal  points  N,  was  varied  through 
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10,  20,  30,  40  and  results  for  stress  intensity  factors  are  compared 
in  Table  II  for  a/b  =  .2;  .6.  A  rapid  convergence  may  be  observed 
for  n  >10.  Using  results  for  40  nodes  as  a  reference,  10  nodal  points 
give  about  63%,  and  20  nodal  points  give  94%  of  the  reference  value 
for  a/b  =  .2.  For  a/b  =  .6  the  results  for  n  =  20  is  98%  of  the 
result  obtained  using  n  =  30 .  The  convergence  for  longer  cracks  seems 
to  be  better,  since  the  relative  error  of  load  fraction  distributed 
over  the  crack  (governed  by  element  size)  is  smaller  for  long  cracks. 

The  stress  intensity  factor  solution  obtained  for  the  uniform 
discretized  stress  loading  on  the  central  segment  is  given  in  Table 
III,  and  is  compared  to  published  results  for  a  uniform  stress  loading 
(Ref  20:11).  Table  III  shows  that  present  results  are  low  by  five 
percent,  a  factor  that  may  be  attributed  to  the  number  of  nodal  points. 

In  order  to  check  the  overall  approach. cases  of  tc/to  =  1  were 
run.  This  case  is  actually  a  case  of  a  center  crack  gage  with  uniform 
thickness  and  half  length  of  Lc  +  L0,  subjected  to  prescribed  displace¬ 
ments  at  the  end.  Comparison  to  published  results  (Ref  20:16)  indicated 
that  present  results  are  within  10-15%  higher  than  published  for  the 
same  aspect,  ratios,  and  uniform  thickness. 

Due  to  the  fact  that  the  gage  is  constructed  from  relatively  short 
plates,  the  finite  elements  solution  was  checked  thoroughly,  and  re¬ 
sults  were  found  reasonable  (see  Appendix  C) .  Hcwever,  the  finite  ele¬ 
ment  analysis  was  performed  on  a  complete  outer  segment  instead  of  a 
half  since  the  displacements  results  were  more  reasonable  in  this  con¬ 
figuration. 

During  the  center  segment  analysis,  attempts  to  invert  the  flexi¬ 
bility  matrix  were  made  and  difficulties  encountered.  However,  no 
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Table  I 


Convergence  of  Stress  Intensity  Factor 


L 


c  ~ 


.4 


b  =  2.0 


n  =•  20 


M**  -  12 


M  =  18 


M  =  24 


1.2620 


1.2660 


1.2662 


*n  =  number  of  nodal  points  on  half  width 

**M  =  number  of  terms  in  expansion 

2M  =  number  of  coefficients  in  expansion 


Table  II 


Convergence  of  Stress  Intensity  Factor 

k x/osA> 


L 


c 


.4 


b  =  2.0 


M  =  18 


a/b  n  =  10  n  =  20  n  =  30  n  =  40 

.2  .8497  1.266  1.322  1.30 


o 


6.607 


6.731 


t 

I 

1  Table  III 

1  Stress  Intensity  Factor  for  Stress  Loading  -  Comparison 

with  Published  Results 


* 

t  ' ; 

r  i 

a/b 

Lc/b  =  .5 

Kj/OgA) 

Present  Results 
Discrete  load 

Rcoke  (Ref  20) 
(Data  transferred  to 
same  nondimensionai 
form) 

.2 

.84 

.88 

* 

.4 

1.70 

1.78 

.6 

3.15 

3.28 

explanation  for  this  difficulty  was  found.  The  final  results,  as  well 
as  the  intermediate  results  of  the  analysis,  were  within  reasonable 
agreement  with  other  published  results.  Thus,  the  flexibility  matrix 
was  accepted  in  spite  of  the  unexplained  difficulty  in  computing  an 
inverse . 

Figures  9  through  19  show  tile  stress  intensity  factor  results 
obtained  in  this  study.  One  of  the  gages  (Lc/b  =  .2  Lo/b  =  .63 
tc/tQ  =  .5)  was  developed  by  McDonnell  Douglas  Corporation  as  part 
of  the  evaluation  of  the  crack  gage  on  F-4  full  scale  fatigue  test 
(Ref  14)  .  Another  gage  (Lc/b  -  .625  L^-Zb  =  .625  tc/to  -  .2)  was 
developed  by  Boeing  Wichita  as  part  of  contract  F33615-77-C-507 3  (Ref 
13) .  The  results  are  given  in  a  nondimensionalized  form  (Kj/Og/b) 
(Es/Eg)  vs.  a/b,  where  as  is  the  stress  of  tie  carrying  structure 
at  the  bond  line,  such  that  a  longitudinal  displacement 

U*  =  — -  Lt  (43) 

E  1 


\ 

.1 
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is  transferred  to  the  gage ,  b  is  the  width  of  the  gage  and  a  is 
the  length  of  the  crack.  Figure  9  shows  the  typical  benefits  of  the 
stepped  gage  compared  to  a  gage  with  uniform  thickness.  For  the  same 
aspect  ratio,  i.e.,  1.5,  the  stepped  gage  gives  a  40%  higher  stress 
intensity  factor  and  is  more  constant.  Furthermore,  the  range  of 
stress  intensity  factor  available  to  choose  for  a  specific  need  with 
geometrical  restrictions  is  very  large.  Figures  10  through  18  show 
that  stress  intensity  factor  will  grow  as  thickness  ratio  will  de¬ 
crease,  and  as  L0  grows  for  a  given  Lc  .  Figure  19  is  typical, 
showing  the  effect  of  change  in  thickness  ratio  change  while  the 
length  ratio  is  held  fixed.  It  may  be  observed  that  changes  in  length 
ratios  of  center  and  outer  segment  appear  to  change  the  magnitude  and 
the  shape  of  S IF  vs.  a/b,  while  thickness  changes  will  affect  primarily 
the  magnitude  only.  It  was  found  that  for  a  thickness  ratio  less  than 
tc/tQ  =  .1,  the  stress  intensity  solution  is  as  if  the  prescribed  dis¬ 
placement  was  applied  directly  to  the  center  segment. 

For  fast  approximations  of  the  stress  intensity  factor,  the  following 
equation  may  be  used  (see  Appendix  D  for  details) . 


Lo  +  Lc 


\/b 


Tf  b 


£d.  ie 


/E-a*  +  /(1-v2)L 


(44) 


In  this  equation,  Rice's  (Ref  22)  limit  values  for  Kj  in  a 
short  infinite  strip  with  a  semi-infinite  crack  is  encountered.  The 
values  of  Kj  obtained  by  this  approximation  are  within  10-15%  lower 
than  the  actual  solution,  for  a/b  >  .6.  The  approximation  was  verified 
for  values  of  tc/tQ  and  L^/Lq  in  the  range  investigated,  and  is  suggested 
as  a  rapid  means  of  obtaining  a  trial  design  with  desired  stress  intensity 
factor. 
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Kr/as/b)  •  (Eg/E 


(Kj/asvi))  •  (Es/E  )  -  Stress  Xntensi 


(Kj/Og/b)  •  (Eg/E 


/b) * (Es/Eg)  -  Stress 


V.  Experimental  Results 

The  methodology  for  the  verification  of  analytical  results  of 
stress  intensity  factor  for  center  cracked  stepped  gage  is  based  on 
Paris  '  law . 

da/dN  =  C(AX)m  (45) 

Where  a  is  the  crack  length;  N  the  number  of  loading  cycles;  C,  m 
the  material  constants  and  AK  is  the  difference  of  stress  intensity 
factor  for  each  loading  cycle.  If  we  measure  da/dN  in  an  experiment, 
we  can  find  AK  ,  provided  the  material  constants  C  and  m  are 
known.  When  conducting  a  constant  amplitude  stress  cyclic  test  with 
R=0  where  R  is  Cmin/csmax  then  AK  becomes  actually  the  stress  inten¬ 
sity  factor.  This  sort  of  test  is  to  be  designed  such  that  no  retar¬ 
dation  effects  will  be  experienced  during  the  test.  A  series  cf  experi¬ 
ments  were  designed  to  check  the  analytical  results  of  this  study  and 
in  addition  were  a  part  of  the  evaluation  of  the  crack  gage  designed 
by  McDonnel  Douglas  for  the  F-4  full  scale  fatigue  test.  Due  to 
technical  difficulties,  this  test  program  is  not  completed  and  inter¬ 
mediate  results  are  not  reported.  However,  Boeing  (Wichita)  as  part 
of  contract  No.  F33615-77-C-5073  (Ref  13)  has  conducted  a  series  of 
base  line  tests  and  constant  amplitude  stress  crack  growth  test  for 
a  similar  stepped  crack  gage.  Since  their  data  was  applicable  for 
verification  of  the  results  presented  in  this  study,  raw  data  was  ob¬ 
tained  and  reduced  to  compare  to  analytical  results  to  be  presented. 

The  material  used  for  the  Boeing  specimens  was  7075-T651  aluminum. 

-14 

The  Paris'  law  constants  were  found  to  be  C  =  .1479*10  and  m  =  2.49145 
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where  the  units  of  da/dN  and  Ak  are  in/cycle  and  psi*Jin,  respectively, 
base  line  data  was  fit  to  equation  (45)  with  a  computer  program  pro¬ 
vided  by  Dr.  Ashbaugh  and  Dr.  Grandt  from  AFML  (Ref  23) .  This  program 
uses  Feddersen's  formula  modified  by  Tada  (Ref  21:2.2)  to  evaluate 
K.j  for  a  center  crack  coupon.  The  gage  shown  on  Figure  20,  was  bonded 
to  the  carrying  specimen,  shown  in  Figure  21,  by  FM-73  adhesive.  Con¬ 
stant  load  amplitude  tests  were  conducted  with  (Js  =  10,000  psi  and 
R  =  0.  Crack  growth  vs.  cycles  was  recorded.  From  a  vs.  N  data,  da/dN 
was  calculated  and  by  applying  Paris'  law,  Kj  was  determined.  The 
results  shown  in  Figure  22  were  taken  from  two  different  tests  employing 
six  gages ,  and  where  gages  were  located  back  to  back  on  the  carrying 
specimen.  The  first  test  carrying  specimen,  AFCG-1,  included  two  of 
the  gages.  On  the  second  test  specimen,  AFCG05 ,  four  gages  were  bonded 
to  the  carrying  specimen.  The  graphical  results  of  Figure  22  show  very 
good  agreement  between  test  and  prediction.  It  has  to  be  noted  that 
there  is  significant  scatter  in  the  experimental  results,  as  is  to 
be  expected.  Scatter  is  inherent  in  this  experimental  procedure  due 
to  variation  in  materials,  loading  accuracy,  crack  length  reading  and 
differentiating  procedure  used  to  obtain  da/dN. 


(Kj/Og^)  (Es/Eg)  -  Stress  intensity 


VI.  Conclusions  and  Recommendations 


This  study  provided  methods  for  obtaining  stress  intensity  factor 
for  trapezoidal  planform  gage  and  stepped  gages.  For  the  stepped  gage 
solution,  a  method  was  developed  using  Reissner's  principle  along  with 
finite  element  techniques  which  proved  to  give  satisfactory  results, 
as  verified  by  test.  It  was  clearly  demonstrated  that  the  stepped 
gage  can  be  designed  to  have  stress  intensity  factors  within  a  wide 
range . 

Conclusions 

1.  The  trapezoid  shaped  edge  cracked  gage  has  no  significant 
advantage  over  rectangular  gages. 

2.  The  method  of  solution  introduced  for  stepped  gages  has  been 
shown  to  be  suited  to  the  analysis  of  these  geometries  with  good 
accuracy . 

3.  The  stepped  center  crack  gage  provides  a  convenient  moans  for 
obtaining  desired  stress  intensity  factors  due  to  the  relative  ease 
with  which  substantial  changes  in  Kj  may  be  obtained  by  varying  the 
geometry,  and  is  to  be  preferred  to  other  geometric  shapes  which  have 
been  considered  for  the  crack  gage. 

Recommendations 

1.  It  is  recommended  that  the  analytical  results  presented  in 
this  study  be  verified  by  further  testing  of  stepped  gages  of  other 
geometries . 

2.  It  is  recommended  that  the  Kj  solution  method  for  stepped 
gages  introduced  be  included  in  crack  growth  prediction  models  so  as 


to  provide  an  analytical  model  for  crack  growth  prediction  for 
crack  gages  under  load  sequence  more  complicated  than  constant 
amplitude . 


-j 
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Appendix  A 


Stresses  and  Displacements  as  Result  of 
Williams'  Stress  Function 


For  the  symmetric  type  problems  introduced,  only  the  even  part 
of  Williams'  stress  function  is  of  interest  and  may  be  expressed  as : 

Xe  =  <-in~lA2n-lrn+1/2{"  cos(n  "  3/2>0  +  ----- -3  cos  (n  +  1/2)0}  + 

2n  +  1 


+  (-1) "a  rn+1{-  cos (n  -  1) 6  +  cos (n  4  1)0} 
2n 


arr'  a00'  °ro  are  found  by: 

1_  32Xe  1  ^Xe 

rr“r2  W~  r 


a  =  8  Xe 

°00  — f 

9r 

1  ^Xe  1  ^Xe 

°r0  "  r  3r0Q  +  r2  96 


.  n-1  n-3/2 


°rg (l)  -  (-1)  r1  '  A2n-1(n  -  3/2) • (n  -  1/2)* 


•  •  {sin(n  +  1/2)6  -  sin(n  -  3/2)0} 


°rQ  (2)  =  (_1>nrn~lA2n{n){(n  +  l)sin(n  +  1)0  -  (n  -  l)sin(n  -  1)0}  (51) 


°r0  °r0 (1)  +  ar0(2) 
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(n  +  1/2) ] 'cos (n  -  3/2)0 


rr  (1) 


_  .  ,  .n-1  n-3/2  , , . 

-  (_1)  r  A2n-l‘{[(n- 


3/2) 2  - 


-  (n  -  3/2) (n  -  1/2)  cos (n  +  1/2)6} 


(53) 


rr(2) 


=  (-l)nrn  1A2n*{[(n  -  1) 2  -  (n  +  1) ] cos (n  -1)0 


n(n  +  l)cos{n  +  1)0} 


(54) 


°rr  °rr(l)  +  °rr(2) 


(55) 


.  . . n-1  n-3/2  ,n  -  3/2 

O00  -  (“D  r  A2n-l(n  +  V2)(r»  -  l/2)-{— +  1/2  cos  (n  +  1/2)9 


-  cos (n  -  3/2)0} 


(56) 


,  n  n-1 


CT0O  (2)  =  r  A2n(n)  +  1^cos^n  +  1)0  -  cos  (n  -  1)  } 


(57) 


°00  G09(i)+  °ee(2) 


(56) 


Solving  the  equation 


2  9 

V  Xe  "  8r  (r  3G_) 


(59) 


get  as  a  result 


,,  _  ,  ..n-1  r  3/2  ,  r  (n-3/2) 2  -  (n  +  1/2) 2 

'Pe  “  (-D  - •—  A  *1 - sin  (n  -  3/2)0]  + 

n  -  1/2  n  -  3/2 

+  (~Dn  -  •  [•■■■—■ —  "-(n  — —  sin (n  -  1)0]  (60) 

n  ^ n  n-1 


Substitute  in 


2MUr  -  ~ ~  +  (l-s)r^ 
r  dr  30 


(61) 
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(62) 


2V1U0  -  -  7  30  +  (1  •  s)r  37 

v  =  Poisson  ratio  s  =  v/(l  +  v) , 

)J  =  shear  modules. 

As  a  result: 

2yur(i)  =  (-l)n-V-1/2A2n_1{[(n  +  1/2)  +  .[<n  -  3/2) 2  -  (n  -  1/2)1]  ' 

in  -  1/  2) 

•cos (n  -  3/2)6  -  (n  -  3/2)0cos(n  +  1/2)0}  (63) 

2^Ur  (2)  =  ("l)nrnA2n{  [(n  +  1)  +  ~  S)  »  [(n  -  1) 2  -  (n  +  l)2]]* 

•cos  (n  -  1)  e  -  (n  +  l)cos(n  +  1)0}  (64) 

2pUr  -  2u(U^^j  +  )  (65) 

2lJU0(l)  '  (-«r‘1rn‘VJ*2n.1{  Hn  -  3/2)  ♦  ' 

fKn  "  3/2) Z  -  (n  +  1/2) )]  sin(n  -  3/2)0  +  (n  -  3/2)  sin  (n  +  1/2)0} 

(66) 

2yU9(2'  =  ^^^n^  ~(n  “  +  1  n  S  l  (n  -  l)2  -  (n  +  l)2]]' 

•sin  (n  -  1)0  +  (n  +  1)  sin(n  +  1)0}  (67) 


54 


- - — —  - — 1 

1 

Appendix  is 


Description  of  Computer  Program 

Three  computer  programs  used  to  generate  the  stress  intensity 
factor  in  this  study  will  be  described  here.  Program  COEFDT  was  used 
to  solve  the  trapezoid  case.  COEFINA  was  used  to  provide  the  solution 
for  the  center  segment  of  the  stepped  gage.  Program  TAYLOR  was  used 
to  match  the  segments  of  the  stepped  gage.  A  few  subroutines  are 
common  to  both  COEFDT  and  COEFINA..  The  subroutines  and  parts  of  COEFDT 
and  COEFINA  were  developed  by  Torvik  for  a  previous  analysis  (Ref  9, 

10)  . 

Subroutine  DPS I 

In  this  subroutine  tangential  displacements  in  Eqs  (66)  ,  (67)  are 
evaluated.  The  results  are  transferred  to  the  main  program  as  for 
Eq  (66)  and  Eq  (67)  separately. 

Subroutine  DR 

In  this  subroutine  radial  displacements  as  in  Eqs  (63)  ,  (64)  ara 

evaluated. 

Subroutine  SP.R 

This  subroutine  evaluates  radial  stresses  as  in  Eq  (53) ,  (54) . 

Subroutine  STT 

This  subroutine  evaluates  tangential  stresses  as  in  Eqs  (56) ,  (57) . 

Subroutine  SRT 

This  subroutine  evaluates  shear  stress  as  in  Eqs  (50) ,  (51) . 

Program  COEFDT 

This  program  evaluates  the  matrix  Amq  as  in  Eq  (24)  .  The  eval- 
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uation  is  done  for  each  of  the  three  edges  of  the  three  separately, 
and  the  results  added  to  produce  the  required  matrix  element.  The 


evaluation  is  made  first  on  the  uncracked  edge,  and  last  on  cracked 
edge.  Then  Fq  is  evaluated  on  edge  2  via  Eq  (25).  Finally,  the 
system  of  Eq  (23)  is  solved  for  c^,  and  dp  is  the  desired  result. 

Program  COEFIEA 

This  program  evaluates  the  matrix  A^g  via  Eq  (29)  in  the  same 
order  as  done  in  program  COEFDT.  The  vector  Fq  is  evaluated  on  edge 
2  via  Eq  (30)  by  moving  a  unit  load  from  one  nodal  point  to  the  other. 
Fq  is  evaluated  n  times  as  the  number  of  nodal  points .  Solution 
of  Eq  (23)  each  time  gives  finally  an  array  {d}  of  coefficients  dr,  . 
In  parallel  to  the  above  described  procedure,  the  program  constructs 
the  flexibility  matrix  [Sc]  for  the  center  segment. 

Program  TAYLOR 

This  program  evaluates  Kj  via  Eq  (40) . 
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Appendix  C 


Fini. te  Element  Solution 

The  flexibility  matrix  for  the  outer  segment  of  the  gage  was 
obtained  by  using  a  slightly  modified  version  of  the  finite  element 
analysis  program  "ANALYZE"  (Ref  18) .  The  problem  was  solved  with  a 
mesh  that  consisted  of  452  nodes  and  441  elements.  The  elements  used 
were  quadrilateral  and  triangular  membranes .  Mesh  organization  and 
boundary  conditions  are  shown  in  Figure  23.  As  a  check  of  the  solution, 
displacements  of  the  "pulled"  edge  were  examined.  The  displacements  at 
the  edges  looked  extremely  high  in  comparison  with  displacements  at 
the  center,  but  were  found  to  agree  with  a  second  solution  obtained 
through  NASTRAN.  Since  in  our  case  we  have  a  very  short  membrane, 
width  =  4,  L  =1.9;  1.22;  1.6;  2.0;  4.0;  and  tQ  =  .08,  it  was  suspected 
that  the  restraint  in  the  x  direction  has  an  effect  of  singularity. 

To  check  this  point,  the  restraints  in  the  x  direction  at  nodes  1, 

441  were  released  and  then  all  restraints  in  x  direction,  except  at 
the  mid-node  221,  were  released.  This  modification  did  not  change 
the  results  significantly,  as  can  be  seen  in  Figure  24.  This  suggests 
that  no  singularity  effects  were  introduced  by  the  x  direction  restraint. 
Figure  24  also  shows  the  displacement  variation  as  the  length  of  the 
membrane  grows . 
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Finite  Element  Model  and  Mesh  Organization 
for  Outer  Segment  of  Gage 


Appendix  D 


t 

( 

Stepped  Center  Gage  -  Stress  Intensity 
Factor  Approximate  Equation 


( 


Figure  25.  Schematic  Stress  Distribution  at  Crack  Tip 

of  Stepped  Gage 


Using  the  relationship  between  stress  intensity  factor  and  stress, 
on  the  line  y  -  0,  we  get: 


a  - 


/2ttx 


(69) 


and  the  force  Pa  produced  over  the  uncracked  region  will  be 
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Pa  =  v'f  fcc  KI 


(70) 
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This  force  will  produce  an  average  displacement  on  the  outer  segment 
of 


u 
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P  L 
a  o 


A  E 
o  o 


P  L 
a  o 


b*toEo 


(72) 


The  displacement  of  the  center  segment,  may  be  approximated  by  using 
Rice's  (Ref  22:249)  result  for  an  infinite  strip  with  semi-infinite 
crack. 
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(73) 


The  total  displacement  may  be  expressed  as 


tot 


uc  +  u0 


(74) 


On  the  other  hand,  UtQt  produced  by  the  displacement  of  the 
bond  line  with  the  carrying  structure  is 

Vs 


U 


tot 


(75) 


Substituting  and  rearranging  and  making  a  nondimensional  stress 
intensity  factor,  normalized  with  respect  to  the  nominal  stress  in  the 
structure  we  get: 


fl  _  Lo  *  Lc 

0gA>  L0  tc  /b-a'+  /(1-v  2^c) 


(76) 
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